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Abstract. A limit theorem for a sequence of diffusion processes on graphs is proved in a 
case when vary both parameters of the processes (the drift and diffusion coefficients on every 
edge and the asymmetry coefficients in every vertex), and configuration of graphs, where the 
processes are set on. The explicit formulae for the parameters of asymmetry for the vertices 
of the limiting graph are given in the case, when, in the pre-limiting graphs, some groups of 
vertices form knots contracting into a points. 



Introduction 

The main object considered in the present paper is diffusion processes on graphs; a graph is 
treated as a one-dimensional topological space with branching points, rather than as a discrete 
scheme. Such processes arise naturally, on the one hand, in a number of applied models (e.g., 
in a model describing the motion of nutrients in the root system of a plant, see pQ) and, on the 
other hand, in some theoretical constructions (e.g., in the study of small random perturbations 
of Hamiltonian dynamical systems, see [2], or in the study of the asymptotic geometric properties 
of discrete groups, see [3|). Such processes possess a number of interesting and nontrivial internal 
structural peculiarities. Let us mention one of them, that was revealed by B.S.Tsirelson and 
consists in the following (see detailed review in [3]): the typical diffusion on a graph generates 
a filtration, that cannot be obtained from the filtration generated by some (possibly, infinite- 
dimensional) Brownian motion, "in real time", i.e., by means of a morphism. 

The full description of a diffusion process on a graph in terms of its infinitesimal operator is 
given in [5] . In a nonformal way, such a process can be described as a mixture of the motions 
"along an edge" and "in a neighborhood of a vertex". A motion of the first type is described by 
a one-dimensional diffusion process and is defined by its coefficients of drift and diffusion. To 
describe a motion of the second type, it is necessary to set additionally the parameters playing 
the role of boundary conditions at a vertex, that define the behavior of a process in the vertex. 
In [5], these objects are called "the gluing parameters". We call them also "the asymmetry 
parameters", since the construction of the process, that we use as the basic one, differs from 
the one developed in [Sj. The purpose of this paper consists in the description of the limit 
behavior of a diffusion process in the situation where both the above-mentioned parameters of 
the process and the graph itself vary. If the limiting lengths of edges of the graph are nonzero 
(i.e., the vertices of the graph do not "glue together", and the configuration of the graph, in 
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fact, does not vary), then the required limiting result is very similar to the standard limit 
theorems of the theory of diffusion processes. The new specific problems requiring a separate 
analysis arise in the case where some groups of vertices form knots contracting into a points. 
In this case, the description of the limiting behavior of diffusions can be investigated in the 
framework of the approach developed in [5J, and non-trivial problem consists in calculation of 
the asymmetry parameters of the limiting process. As an example of the expressions that can 
be finally obtained, we give the model, in which the graph consists of two vertices and three 
edges: the edge with length e joining the vertices and two half-lines beginning at these vertices. 
If the diffusion on each edge is a Brownian motion, then its distribution converges, as e — * 0, to 
the distribution of a skew Brown motion on a real line ([£]). Its asymmetry parameters are set 
by a single "skewing parameter" q (p± = which is calculated in the given model as follows: 
if qi,q2 are analogous parameters for the vertices of a prelimiting graph, then q = th(Ai + A2), 
where Ai^ = arcthgi^, the = ^ e ~^_ e is a hyperbolic tangent, and arcth is the function inverse 
to th. This simple example shows that the determination of the asymmetry parameters for a 
vertex, that is obtained as a result of the contraction of a group of vertices into a single point, is 
a nontrivial problem. In the present paper, we propose a method of the solution of this problem. 

Let us describe one possible application of the main result of the paper. Earlier, we have 
mentioned the Tsirelson's result that states that if, for some vertex of a graph, at least 3 it's 
asymmetry parameters are nonzero (i.e., the vertex of a graph has multiplicity > 3 and is a triple 
point), then the diffusion X is essentially singular in the following sense: there does not exist 
any morphism of the filtration, generated by any Brownian motion, to the filtration generated 
by X ([4]). We call further the singularity, related to the presence of a triple point (in the 
above-mentioned sense), the Tsirelson's singularity. 

The theory of diffusion processes includes a number of results, in which the diffusions, con- 
taining singularities of some type (such as an asymmetric semipermeable boundary at some 
point, or a partial reflection with delay at the boundary of a domain), are represented as the 
weak limits of nonsingular diffusions (see [7], [8]). In these results, typically, the coefficients of 
prelimiting diffusions, in a certain sense, model the singular terms that are present in the lim- 
iting process (such as the singular drift coefficient a = q5o of a skew Brownian motion). The 
Tsirelson's singularity is related to the structure of the phase space (the presence of a nontrivial 
branching point), rather than to properties of the coefficients of the process. The main theorem 
of the present paper allows us to represent a process, possessing the Tsirelson's singularity, as 
the weak limit of the processes without such singularities. We will construct such a representa- 
tion, changing the phase space of the process and representing a vertex with multiplicity > 3 
as a result of "the contraction into a point" of a knot, whose vertices have multiplicity < 2 (see 
Example 2 below). Such a trick gives, in perspective, the possibility to study the properties of 
stochastic flows (i.e., the processes describing a motion of a families of points) corresponding to 
the diffusion on a graph, that possesses the Tsirelson's singularity and for which, by this reason, 
one cannot define the common law of motion of the family of points as the strong solution to a 
system of SDE's. 



1. Basic notation and constructions 



1.1. Phase space. Everywhere in what follows, a graph means a connected metric space 
composed of a finite number of subsets (edges of the graph) homeomorphic to a segment or a 
half-line. It is assumed that the only intersection points for these subsets are the images of the 
ends of a segment or a half-line; all such points are the vertices of the graph. 

It is convenient to consider the graph to be oriented, by assuming that if a point, whose motion 
is described by a diffusion process, is not positioned on a definite edge, then this point can move 
into this edge, only by passing through the separated vertex of the edge (its "beginning"). 
It is clear that this does not restrict the generality, because any non-oriented edge-" segment" 
can be represented as two oriented copies; for an edge-" ray", its beginning, obviously, is the 
single vertex belonging to it. For the vertex i of a graph 9, we denote, by Tj, the family 
of vertices joined with it by edges beginning at the vertex i, and these edges are denoted as 
Llj,j G "?i,r = 1, ... ,R(i,j). The necessity to introduce the additional parameter r is caused 
by that two vertices can be joined by several edges. Analogously, the edges that leave the vertex 
i and are homeomorphic to a half-line are denoted as L r ioo ,r = 1, . . . ,i?(i,oo). The interior of 

the set UjelPi Ur ^1 j U Ur-^ioo ^ s denoted as Oj and called the maximum neighborhood of 
the vertex i. 

On each edge, we introduce a natural parametrization in the following way: for the edge L r ioo , 
we consider the homeomorphism with [0, +oo) to be fixed and define the coordinate of a point 
on the edge as the number corresponding to it via this homeomorphism. The coordinates of 
points on the edges L\ j,j £ Tj are determined analogously, by considering the homeomorphism 
of the edge with [0, Z£ •] to be fixed. 

We make no assumptions that the graph is imbedded in any larger metric space (e.g., that it is 
planar). On the other hand, the natural parametrization allows us to homeomorphically imbed 
each of the maximum neighborhoods Oj in M 2 , the corresponding image being a part of a bundle 
of half-lines in M?. If this does not cause misunderstanding, we will omit the corresponding 
homeomorphism, by considering that any Oj is a part of a bundle of half-lines in M 2 and that 
any edge is a part of M + . 

1.2. Construction of the process. In our consideration, we will use two constructions of 
a diffusion process on a graph. One construction sets the process in terms of its infinitesimal 
operator and is given in [5]. The other, more explicit construction describes the diffusion process 
in terms of its excursions. We cannot give a reference, where the required construction would 
be described in the generality sufficient for our purposes, therefore, we give its description here. 
Of course, we do not pretend for a priority, because various versions of such a construction 
were given earlier (see [9], §4.2 and |10].|11|). and its main idea is widely known. In the next 
subsection, we will show that both versions of the construction of a diffusion process on a graph 
lead to the same object. 

First, let a graph S be a bundle of half-lines, i.e., it has one vertex and R edges-" rays" . 
The Walsh's Brownian motion is set on such a graph as follows (|10|.|llj): it is a continuous 
Feller Markov process; on the start from a point on one of the rays, up to the first moment 
to hit of the vertex, this process propagates along this ray as the Wiener process. On the 
start from the vertex, the process is constructed as follows: one take the Brownian motion B 
on [0, +oo) with reflection at zero, and construct a sequence of identically distributed random 



values {e n ,n > 1}, that are independent of the process B and one another and such that P{e\ = 
k) = Pk>k = 1, . . . , R, where pi, ■ ■ ■ ,pr are the given numbers ("the asymmetry parameters" of 
the vertex), pi + • • • + pr = 1. Then one enumerate, in an arbitrary measurable way, all the 
excursions {[a m , b m ]} of the process B at zero and denote, by m{t) for an arbitrary time moment 
t, the random variable setting the number of a current excursion, i.e., such natural number that 
a m(t) < t < b m u\. Such a variable is well defined on the set {B t > 0} having probability 1. Now 
the value of the required process X at a time moment t is determined in the following way: it 
is located on a ray with the number e m M, and its coordinate on this ray is equal to Bt. 

Next, let the graph 9 be the same one as above, and let, for each edge U , the functions a r , a r 
on the ray (0, +oo), which are interpreted as the coefficients of drift and diffusion of the process 
on the edge, be set. We assume that these functions are measurable, the functions a r , [c r ] are 
locally bounded, and a r is locally integrable. In addition, in order to shorten the consideration 
and to exclude the possibility of an "explosion", we assume that the functions a r ,a r possess at 
most linear growth at infinity. We introduce a new parametrization 

f x f y 2a r (z) 

x = S r (x) = exp[- / - dz] dy, x G L r , (1) 

Jo Jo F \ z )\ 

and, for the given collection p = {p\ , . . . , pr} of the asymmetry parameters of a vertex, construct 
a Walsh's Brownian motion X with these parameters. We introduce a process {r(s), s > 0} such 
that X s G L r ( s \ s > 0, almost surely and put 

e t = J* {a r ( s \x s )-[s r ^}'(x s )yds, T t = er 1 ^mf{u\e u >t}, t>o. (2) 

Now, we put Xt = S~ 1 (X Tt ),t > 0, where S" 1 is the change of a parametrization on the initial 
graph which is inverse to S. By construction, X is the Feller Markov process with continuous 
trajectories. The proof of this assertion is exactly the same as that for an analogous proposition 
about the general one-dimensional diffusion (see [9], Chap. 3). We interpret process X as the 
diffusion process on a bundle of half-lines with the coefficients of drift and diffusion on rays 
{a r , a r } and the asymmetry parameters of the vertex p. 

At last, let a graph 9 be arbitrary, let the collection of asymmetry parameters pi be set for 
its each vertex, and let the coefficients of drift and diffusion {a^of^} be set for each edge 
L$j. Supposing, e.g., that a r id (x) = a r id (x) = 1, j ^ oo,x > l r i;j (or a r id (x) = a^-j), ^{x) = 
a lj{l<ij)-,3 ^ oo,x > l\ j, if a\^o\j are continuous on [0, /[■]), we can consider that these 
functions are given on the whole M + . Now, the Feller Markov process with continuous trajectories 
is well defined on 9 via the following convention: on the start from a point lying in the maximum 
neighborhood Oj of any vertex i, it moves as a diffusion process on the bundle of rays with the 
parameters pi, {a\ p £ {1, . . . , N} U {oo}, r = 1, . . . , R(i,j)} till the exit from Oj. At the 

moment of exit from Oj, it is located at some vertex j G Jj, and, after this moment and up 
to the next exit time (from Oj), it moves as a diffusion process on the bundle of rays in the 
maximum neighborhood Oj, and so on. The process, constructed in the way described above, 
is the diffusion process on 9 with the parameters {pi,a^ p o\ j}. 

Let us note that the constructed process spends zero time at every point of the phase space 
with probability 1. A wider class of the processes with "sticky" points can be constructed by 



means of a time change (see the details, e.g., in [12] . Section 3.3). In order to shorten the 
exposition, we exclude such processes from consideration. 

1.3. Martingale description and infinitesimal characteristics of the process. 

To describe the infinitesimal characteristics of the process constructed in the previous section, 
it is sufficient to consider the case where the graph S is a bundle of half- lines. For such S, let 
the diffusion X with the parameters p = {p r }, {« r , o~ r } be given. We assume that the functions 
a r ,a r , in addition to the above-imposed conditions, are continuous on [0, +oo). By X r , we 
denote a coordinate process on the r-th edge: X\ is equal to the coordinate X t on the edge L r , 
if Xt lies on this edge at a time moment t, and X[ equals zero otherwise. A vertex of the graph 
is denoted by the letter O. The process X can be described in terms analogous to those of the 
Skorokhod problem for a process on a half-line with reflection at a point of the boundary. 

Proposition 1. There exists a nondecreasing process Vt such that 

(i) it increases only when the vertex is visited by the process X, i.e., 

I l{X s ^O} dV s = a.s., t> 0; 
J o 

(ii) for an arbitrary edge L r , the process 

Ml = Xl-p r V t - f a r {Xl)l [x ^ Lr} ds 
Jo 

is a continuous martingale with the quadratic variation 



{M r 



>t 



[\a r (X r s )fl {XseL r } ds. 
Jo 



Proof. For a Walsh's Brownian motion, this assertion follows from the reasoning analogous to 
that given in [4], Section 3: each of the processes X r belongs to the class S + ([13J, VI. 4. 4.), i.e., 
can be represented as a sum M r + V T , where M r is a local martingale, and V r is a nondecreasing 
process such that J Q lix r >o\dVg = 0. It is easy to verify that the quadratic characteristic of 
M r is equal to J" * 1{xj>o} ds. On the other hand (see [13], VI. 4.4.), 
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Let Vt be the local time of the process X at a vertex O, i.e., 

1 

e'-^o+ 2e j 

It is easy to see that both V r and V are VF-functionals of the process X. Moreover, their 
characteristics f r and / are connected by the relation f r = p r ■ f by construction. Using 
Theorem 6.3 in [T4], we get V r = p r ■ V . That is, Proposition 1 is valid for a Walsh's Brownian 
motion, and the process V is the local time of the process X at zero. 

Now let Xt = S' _1 (V Tt ), where X is a Walsh's Brownian motion (we use the construction of 
the previous section), and Vt is its local time at the vertex. Then Proposition 1 is valid with 
Vt = V Tt , that follows from the ltd formula (note that [5* r ]'(0) = 1). 



Assertions (i),(ii) can be naturally interpreted as a version of the martingale problem for the 
pair of processes (X, V). Another version of the martingale problem for the process A is given 
below. 

Proposition 2. Let a continuous function (ft on S be equal, on every edge L r , to some function 
4> r E C 2 ([0,+oo). Denote A<j>{x) = a r (x)[<f> r ]' (x) + ±[a r (x)] 2 [cj) r ]" (x), x £ U and A (<t>) = 

Then, for an arbitrary cf) satisfying the above-indicated condition and such that Ao(<j>) = 0, 
the process 

Mf = <f>{X t )- [ A<j,(X s )ds 
Jo 

is a continuous martingale. 

This assertion follows immediately via the Ito formula from Proposition 1 and the fact that 
the process X spends zero time at O. 

It is important that, as the following theorem shows, both the martingale problems given in 
Propositions 1,2 are well-posed. By Da, denote the set of continuous bounded functions 4> on 
S, such that the function Acf) is well defined, continuous, and bounded on S, and the condition 
AoO) = holds. 

Theorem 1.1) The operator (A, Da) is a n infinitesimal operator of the process X constructed 
in the previous subsection. 

2) The process X is the unique solution of the martingale problem posed in Proposition 2, 
endowed by the given initial distribution P{Xq £ •). 

3) Let V be the process constructed in the proof of Proposition 1. Then the pair of processes 
(A, V) is the unique solution of the martingale problem posed in Proposition 1, endowed by the 
given initial distribution P(Xq £ •), that satisfies the following conditions: 

(i) P(V = 0) = I; 

(ii) with probability 1, the process X spends zero time at the vertex O. 

Proof. The fact that (A, T>a) is a pre-generator of the process A follows from Proposition 2. 
Theorem 3.1 in [5] ensures the fact that Da is the whole domain of definition of the generator 
of the process A. The second assertion follows from Theorems 3.1 and 2.2 in [S]. The third 
assertion is a consequence of the second one and the Ito formula. 

According to Theorem 1, the constructive description of a diffusion process on a graph, 
presented in the previous section, and the semigroup description presented in Section 3 in [5] 
are equivalent. 

1.4. Limit theorem for a graph with constant configuration. The construction of 
Section 1.2 yields directly the following limiting result for the sequence of diffusion processes 
{X n } on a graph 9 with the parameters {p", a^'J, o~l'™}. 

Theorem 2. Let 

r.n r r /ttt> + \ r i n r 

a i,j ~^ a i,V L hloc{^ ), <Tij -^O-ij 

locally uniformly on [0, +oo), and let 

Pi = {P2j,j G Vi U {oo}, r = 1, . . . R(i,j)} -^ Pl = {p r ltj ,j € 7i U {oo}, r = 1, . . . R(i,j)} 



componentwise (recall that we suppose that the coefficients a r ^,a r ^ are given on [0, +00), by 
setting them, if necessary, to a constant on (l^,+oo)). 

Then the sequence of the distributions of processes {X n } converges weakly in C(M + ,S) to the 
distribution of a diffusion process X with the parameters {pi, a\ •, a\ •}. 

Proof. In each of the maximum neighborhoods Oj, we consider the changes of the phase 
variable S n , [5 n ] _1 , that are defined by the coefficients X n on the corresponding edges. By 
virtue of the imposed conditions, S n , [S" 1 ] -1 converge locally uniformly to S, [5*] _1 , defined by 
the coefficients of X. By performing the change of a phase variable for X n in each of the 
neighborhoods Oj, that is inverse to (1), and the change of the time that is inverse to (2), 
we obtain the process X n being a composition of the Walsh's Brownian motions, switching 
themselves at the time moment of the transition from one neighborhood to another one. The 
integrands {cr n [5 ri ] / } 2 in the time change (2) time also converge on each edge locally uniformly 
to {crfS 1 ]'} 2 . The probability of the event, that at least M transitions from one neighborhood to 
another one occur for the process X n , can be estimated uniformly in both n and the starting 
point by a term of the form Ca M , where a G (0,1), C are some constants (this follows from 

the estimate P(sup t<£ |W(i)| > c) < e e for the Wiener process W; see also an analogous 
estimate for the ltd processes in [15) . Lemma 8.5). Thus, by virtue of the strong Markov property 
for X n ,X, the proof of Theorem 2 is reduced to the following. There exists the sequence of 
the Walsh's Brownian motions {Z n } on a given bundle of rays starting from the vertex, and 
the corresponding collections p n converge componentwise to the collection p for the Walsh's 
Brownian motion Z. There is also a family of neighborhoods of the vertex 0™, lim„ 0" = 0, and 
Tq2 and Tq are the moments of the exit of the processes from n and 0, respectively. We have 
to show that the distributions of the pair (Z n , Tq2) converge weakly to the distribution [Z, Tq). 
This can be proved, by leaning on three following simple assertions. First, Z n converges weakly 
to Z (here, we can explicitly write the transition probabilities for Z n ,Z). Secondly, Tq, for an 
arbitrary neighborhood 0, is an almost surely continuous functional of the trajectory Z. At 
last, Tq is a monotonous and, for almost all trajectories Z, continuous function of 0. Two last 
facts are a consequence of analogous assertions for the Brownian motion. Theorem 2 is proved. 

Remarks. 1. The assertion of Theorem 2 could be proved, by using the general limit theorem 
4.1 in [5j. However, the reasoning presented in the proof is important for us by themselves, 
because it is, in essence, a part of the proof of the main result of the present paper, Theorem 3. 

2. The assertion of Theorem 2 remains valid if we assume that, for each {X™}, its diffusion 
coefficients are written w.r.t. its own parametrization \£ n , and \P" o — > id, ^ o — > id 
locally in C 2 (M + ) on each edge. Informally, this means that the lengths of edges can vary, not 
tending to zero. 

2. Limit theorem for graphs with variable configuration 

2.1. Statement. Further, we assume that the diffusion processes X n are set on the graphs S™ 
with identical combinatorial configuration (i.e., with the identical procedure to join edges), but 
with different metrics. Formally, this corresponds to the setting of different parameterizations 



ty n on the same graph 9- Nonformally, this means that the family of edges is not changed, 
but their lengths vary. We assume that the characteristics {pf, a^, of'j 1 } of the processes X n 
(relative to the corresponding parametrizations Vl/ n ) satisfy the conditions of Theorem 2 and 
consider the question about the limiting behavior of the distributions of the processes X n . 
We assume that the lengths of some edges tend to zero. For i, j such that 3r : ip" — > 
7^ 0, we write i -~ » j. For i,j such that 3ji, . . . , j m : i — ■* — >■ j 2 , • • • , jm — » J, we 
write i j. We require that the following symmetry condition be satisfied: 

i -w j 44> j -w i. (3) 

If condition (3) fails then the limiting process can fail to be a diffusion. Namely, the limiting 
process (in the sense of the convergence of finite-dimensional distributions) can possess discon- 
tinuous trajectories if condition (3) does not hold. 

By construction, "W is the equivalence relation; the collection of vertices connected by this 
relation will be called "a knot" . It is natural to define a new graph 9, on which a limiting process 
will be defined at last, as a graph, in which the vertices are the knots of the initial graph, and 
the edges are those edges of the initial graph that have not contracted into a point. We impose 
the following natural condition: 

A. On each of the edges Lp™ not contracting into a point, the parameterizations Vl/ n converge 
in the sense of the convergence of the C 2 -diffeomorphisms of M + to a certain parametrization 
ty, and ij'" — > V: i > as n — > oo. 

The phase space of the limiting process is the graph 9 with the parametrization ^. To 
formulate the limiting result, it is necessary to set the projection {X n } of the initial sequence 
{X n } on this space (see [5], Section 4). This can be performed in the following way: if X" lies 
on a non-contracting edge of the graph 9, then X" lies on the corresponding edge of the graph 
9, and its coordinate on this edge is obtained from the coordinate of X™ by the transformation 
<]/ o [vl/™] -1 . If X" lies on a contracting edge, whose ends belong to the knot i, then X™ = i. It 
is clear that the trajectories of the process X n are continuous. 

The asymmetry parameters of the limiting process will be determined by the internal struc- 
ture of prelimiting knots, let us introduce the necessary objects and assumptions. The main 
assumption consists in that each knot is homogeneous, i.e., the lengths of all internal edges in 
it tend to zero with the same rate, 

Vi 3 :<^(n)- and Vz,jGi,r = 1, . . . , R(i, j) -f- - V ■ > 0, n - +oo, (4) 

<Pl(n) u 

where the numbers {11 j} are given. This assumption can be weakened, but it cannot be removed 
at all. To shorten the consideration, we omit the details here. We note only that if condition 
(4) is not imposed, we can faced with situations where the limiting process spends a positive 
time at a knot with a nonzero probability. 

We recall that we assume that, for each vertex i, the collections of the parameters p™ converge 
to a certain collection pi = {plj,j € 7% U {oo},r = 1, . . . R(i,j)}. 



Let the knot i be fixed. We set i\A = #{i e i}, 



r<R{i,j) l ' J jev.ni 
Consider the N~ x ^-matrix ^4* defined by 

7, 7, J 



4* = %> M'ei. 

^7 

The matrix A 1 is the matrix of transition probabiiities for some Markov chain. By virtue of 
condition (3), all the states of this chain form unique class of essential states. Therefore, there 
exists the unique invariant distribution for the chain. We denote this distribution by -it 1 and set 
the collection in the following way: each edge U. , of the graph S is represented by some edge 

1,3 

L\- of the graph S with i € i,j i (it is possible that j = oo, then j = oo). For this edge, we 
put 

i>';, ■!>■,■ ( 5 ) 

where P % is the normalizing factor which is defined by the condition 

£ 

ie^U{oo},r<R(7j) 



f., = l. 



At last, we set the functions a~ ~,a~~ as the limits of the functions a r -'"',a^"' on each edge L~~ 

i,3 1,3 l 'J l >3 1,3 

corresponding to the edge L\y Now we can formulate the main result of the present paper. 

Theorem 3. Let the characteristics {p", a['J, of the processes X n satisfy the conditions 
of Theorem 2, and let conditions (3) and (4) be satisfied. Let also the sequence of distributions 
fx n (-) = P(X n (0) £ •) converge weakly to some measure Then the sequence of distributions of 
the processes {X n } in C(M + ,S) converges weakly to the distribution of the diffusion process X 
with the above-set parameters {^,0^,^-.} and with the initial distribution fi. 

2.2. Proof. The reasoning analogous to those used in the proof of Theorem 2 allows us 
to restrict our consideration to the case where a T /l = 0, 0\ r ' n = 1, all the vertices of the initial 

L ,3 l ,3 

graph S form a single knot, and only nontrivial (i.e., not contracting into a point) edges are the 
edges-" rays" that are homeomorphic to half- line. 

Remark 3. For such a reduction, it is significant that the conditions of convergence of the 
coefficients ap™, a^J, which were formulated in Theorem 2, hold also for edges contracting into a 
point. Otherwise, the assertion of Theorem 3 can be violated. For example, if for some (not all) 
edges, contracting into a point, a-'J are constant functions tending to +oo as n — > +oo, then the 
homogeneity condition (4) will be broken after a change of the phase variable. If for some (not 
all) edges, contracting into a point, are constant functions, which tend sufficiently rapidly 
to as n — ► +00, then the time spent by the process X n at the vertex i, will not tend to zero. 

Let us proceed with the proof of Theorem 3 in the above-indicated case. First of all, we note 
that the sequence of distributions of the processes {X n } is weakly compact in C(M + ,9)- The 



simplest way to prove this, is to use the criterion for weak compactness given in [16j . Theorem 
8.2. For an arbitrary e > 0, the process {^t n } outside the vertex neighborhood B(0,e) with 
radius e is the Brownian motion. This easily implies that, on every finite time interval [0, T], 
for the continuity modulus wt(X u ,S) = sup st<T u s \ <s dist(X^,Xg), the following convergence 
holds true: 

sup P{w T (X n ,S) >2s)^0, 5^ 0. 

n 

This allows one to apply the above-mentioned theorem. 

Our aim is to show that any limiting point of a sequence of the distributions of the processes 
{X n } gives a solution of the martingale problem for a Walsh's Brownian motion, that was 
formulated in Proposition 2. Since, by assertion 2) of Theorem 1, this problem is correctly posed, 
this yields the assertion of the main theorem. To prove the required martingale characterization 
of the limiting point, we will study the limiting behavior of resolvents of the processes {-^ n } 
(or, more exactly, the Laplace transforms of their distributions; note that each of the processes 
X n is not Markov) in detail. 

Let the vertex i, the edge-"ray" with the number r < R(i, oo) which leaves this vertex, and 
the function <p G C&([0, +oo)) be fixed. Consider the quantities 

E?(t) = E(cf>(X?),X? G LfJX" = i) = E(<f>(jt?),X? G LFJX " =i), t> 0. 

By Tj, denote the moment of the first exit of the process X n from the neighborhood Oj. For 
Ef'(-), an analog of the renewal equation, written at the moment n, looks as 

E?(t) = Q?(t) + Yl [ t EUt-s)P(X^ = k,r i eds), (6) 
kea>, Jo 

where Qf(t) = E(4>{X?),X? G LT^n > t\XQ = i) (it is clear that Qf(t) > only if i = i). 
Considering (6) for all i, we get the convolutional equation for the vector E n (t) composed of the 
components of E™(t). Let us introduce the Laplace transformations 

/'OO /"OO /'OO 

V?(\)= e- xt E?(t)dt,V?(\)= e- xt Q^{t)dt,Cl k {\)= e^P^ = k,n G dt), 
Jo Jo Jo 

and let U n (X), V n (\), and C n {X) be, respectively, two vectors and a matrix composed of the 
components of U™(\), V i n (X), and C" fe (A). Equation (6) yields 

U n (\) = V n (X) + C n (X)U n (X), U n {\) = [I — C n (X)]- V n (A). (7) 

To describe the limiting behavior of U n , we need the following lemma allowing us to write 
V n (X), C n (X) explicitly. Denote ^(A) = e~ xt <fi(Bt) dt, where B is the Brownian motion with 
reflection on [0, +oo), starting from zero. 

Lemma 1. Let a Walsh's Brownian motion Z with M rays and the asymmetry parameters 
pi, . . . ,pm be given. Let the points z\, . . . , z m be marked on the rays L\, . . . , L m (m < M) at 
the distances l\, . . . , l m from the vertex, and let r be the first moment when Z hits one of these 
points. 



Then, on the start of the process Z from the vertex O, 



/ e- xt E{<p{Z t ), Z t eL k ,r> t) dt = p fc $(A){^pj cth[ljV2\] + Vj] , A > 0, k > m, 

J ° j=l j=m+l 

(8) 

/ e- xt P(Z T = Zk ,T£dt)= Pk cth[l 3 V2X}+ V pj) , A>0,/c<m, (9) 

^0 s/i[^V2A] ^ J 



ci/ic = 4r> c/lc = eC+ 9 6 C , she 



e — e 



She' 2 >^ 2 • 

Proof. It is sufficient to consider the case where, for some 7 > 0, 4>{u) = (f)(0), u G [0,7] (the 
general case can be obtained from it by approximation of general (j) 6 Cft([0, +00) by functions, 
that are constant in some neighborhood of 0). Let B be the Brownian motion with reflection, 
whose excursions have been used in the construction of Z (see Section 1.2). For x > 0, we 
denote, by r x , the time moment of the first passage of the level x by the process B. Then, for 
x < min(7, l±, . . . ,l m ),k > m, 

E(4>(Z t ),Z t G L k ,r > t) = (f>(0)P(Z t G L k ,r x > t)+ 

+ Z~2 / E j)X ((j)(Zt- s ),Zt-s € L k ,T > t - s)P(Z Tx G L 3 ,T X G ds), 
3=1 Jo 

where Ej :X (-) means the averaging over the distribution of the process Z on the start from the 
point located on Lj at the distance x from the vertex. We set P(t x < t) = T x (t). By the 
construction of the process Z, we have 

P(Z t G L k ,r x > t) =p r [l - T x {t)}, P{Z Tx G Lj,T x G ds) = Pj T x {ds). 

For j / k, 

E j)X {<j){Z t ),Z t eL k ,T>t)= [ E(())(Z t - s ),Zt- s e L k ,T>t-s)Q jiX (ds), 

Jo 

Qj jX (s) = P(Wq = 0,0j < s), where W x is the Brownian motion starting from the point x, and 
0j is the moment of its exit from the interval (0, lj) (lj = +cc for j > to). At last, 

E Kx (4>(Zt),Z t €L k ,r>t)= f E(4>(Zt),Z t -s G L k , r > t - s)Q j: x (ds) + F x (t), 

Jo 

where F x (t) = E((/)(Wf), > t), 9 is the moment of exit of W x from the interval (0, +00). Thus, 
we have the convolutional equation for the function H k (t) = E(cf)(Z t ), Z t G L k , r > t): 

H k (t) = p k U(0)[l - T x (t)} + J F x (t-s)T x (ds)\ + Y,Pi J q J H r (t-s-u)Q j:X (du)T x (ds). 

(10) 

We now introduce the Laplace transformations 



rco poo 

G k (X) = / e~ xt H k (t)dt, S X (X) = / e~ M T x (dt), 
Jo Jo 

/*00 POO 

Rj,xW = / e~ xt Qj, x (dt), tf x (A) = / e~ xt F x (t)dt 
Jo Jo 



and rewrite (10) in the following form: 

G fc (A) = ?^[l-5 a (A)]+p A 5»(A)* x (A) + ^p J -5»(A)i? J> (A)G fc (A). (11) 

3=1 

Let us differentiate (11) with respect to x at the point x = 0. Taking into account that (see [9], 
§1-7) 

f S h[(; 7 -x)V2A] 

5 x (A) = i_ i?, x .(A) = <^ shfcV^A] — , 

ch[xV2A] I exp[-xV2X], j>m 



we have 



and, hence, 



[5,(A)]Uo = 0, [^(A)]Uo = <| ^feV2A]' ^- 

-v2A, j > m 



M 



p fc [* :r (A)]U -G fe (A)V2A{^p J cth^V2A]+ Pj}=0. (12) 

j=l j=m+l 

Repeating literally the performed calculations for the process B (which can be interpreted now 
as a Walsh's Brownian motion on a graph with only one edge), we arrive at the equality 

[*MA)]Uo-*(A)V2A = 0. (12') 
Relations (12) and (12') yield (8). Analogously, for H^(dt) = P(Z T = z k ,T £ dt),k < m, we get 

ft M ~t pt—s 

Hl(dt) = Pk / Q\ Jdt - s)T x (ds) + J2Pj / Hl(dt - s - u)Q hX (du)T x {ds), 
Jo j =1 Jo Jo 

where Q\ x {dt) = P(Wg = Ik, Ok £ dt). This equation in terms of the Laplace transform 
G l k (X) = J^e' xt H}{dt) takes the form 

M 

G}.(\) = Pk S x (X)Rl, x (X) + Y^PjS.WRjA^GlW, (13) 

3=1 

where Rl x (X) = e ^ Xt Q\ x (dt) = g^ X ^^| • By differentiating (13) with respect to x at the 
point x = 0, we get equality (9). Lemma 1 is proved. 

We can suppose that 0(n) = i in (4). Since chx = 1 + o(x),shx = x + o(x 2 ),x — > 0, relation 
(8) yields (for a fixed A > 0) 

( r ) 1 

V{ l (\) = <Z>(\)V2\^\ J2 tH +°(-)' ri^+oo, V?(\) = 0, i?L (14) 
In an analogous way, relation (9) yields C n (A) = A n — B n (X), where 

r,n r 
P ',fc Pi ' k r,* - 

T, rT~r~ ^ IV, ^ Pi,' 



r<R(i,k) V2Ar< J R(i,fc) r<fi(i,oo) 1 

4"fc = " B? k (\) = - ; , +o(-), n^+oo. (15) 

^ l r h I 2^he9i,r<R(i,h) IT. ( 



The matrices A n converge to the matrix A introduced before the formulation of Theorem 
3. This matrix is the transition probability matrix for a Markov chain that is ergodic by 
virtue of condition (3). Hence, the matrices A n possess the same property, beginning from 
some n. This allows us to rewrite (7) in the form more convenient for the asymptotic analysis, 
by performing a suitable change of the basis. Namely, we write the decomposition R N = 
(ex) + (tt)- 1 , [R n ]* = (ei)- 1 + (tt), where N is the number of vertices in the initial graph (knot), 
R N and [R^]* are, respectively, the spaces of column-vectors and row-vectors of dimension 
N, e\ = (1, . . . , 1) T , and ir = (wi, . . . , 7tn) is the invariant distribution corresponding to A. 
Next, we choose the basis {e%, . . . , ejv} in (tt) 1 - and the basis {62, • • • 5 &tv} in (^i)~ L m such 
a way that (bi,ej) = 5ij,i,j = 1,...,N, where 5ij is the Kronecker delta, b\ = tt, and, for 
y = (yi,...,y N ) G [R N ]* and x = (xi,...,x N ) G R N , (y,x) = J2k=i x kVk- Writing now the 
matrix (I — A)i & = — (6j, Aj &e&) (i.e., we write the matrix I — A w.r.t. the bases {bi}, {ei}), 

we get a block matrix of the form I — A = ^ , where = (I — A)^, i,k = 2, . . . , N. 

The ergodicity of the matrix A yields that 

[y(I-A) = 0, ye(e 1 ) ± ]^y = 0, [(I - A)x = 0, xe(i)V^0, 

that means that the matrix D is invertible. 

We perform an analogous operation for all n, by introducing the bases {ef},{6™} in such a 
way that (6™,e™) = 5ij, e± = (1, . . . , 1) T , and tt = (tti, . . . , 7Tjv) ? &" = ^ n being the invariant 
distribution for A n . It is clear that this can be performed in such a way that ef — > e^, 6" — > 
— > +00. Then, writing the matrix I — A n w.r.t. the bases {bf} , {e™} , we get the matrix 
'0 
,0 D" 

VP(X) = (b™,V n (\)),U- l (\) = (6™,C/ n (A)) in (7), we obtain 

*"W = ®o ^w)" 1 ^- (16) 

where %(A) = (6», B_"(A)<£>, B»(A) = (%(A))iJ_„ B;'.(A) = (B^fA), . . . , Bijv(A)), 
B n i(ty = (B% i(A), . . . , i?/v j i(A)) T . To invert the block matrix in (16), we will use a suitable 
modification of the Gauss method. We note that, at sufficiently great n, the matrix D n + B n .(\) 
is invertible. Next, by virtue of (15), B^ k (X) = O(^) and 



/ — A n = ( " ^ n ), where D n — * D, [D n ] 1 —> D 1 . Making now the change of variables 



TV 

nS^(A) -» v 7 ^^ 



r 

i,oc 



r<R(i, oo) 
heVi,r<R(i,h) V ih 



> 0, n — > +oo. 



Therefore, the number 7„(A) = Bf x (X) — B™ .(A^-D* 1 + B n . (A)] 1 B™i(A) is nonzero for great n. 
Then, denoting V n {\) = (V?(\), . '. . , V N (X))\ U n {\) = (u%(\), . . .', U N {X)) T , we get (see [IT], 
Lemma 2.3, and, in particular, formula (12)) that, for such n, 



um = bn(A)]- 1 [y?(\) - (Bl(X), [D n + Bl{\)]- l V n {\)) 
U n (X) = [D n + 5 n .(A)]- 1 (-C/ 1 n (A)5 n 1 (A) + V n (X)). 



These relations and the relations 



Vf (A) = O(i), nVi(X) -+ $(A) V2X • 7^ < 



y 1 } , n^+cc, 

i fceT i ,r<i?(i,fc) i> fe 



that are valid by virtue of (14), imply that, as n — > +00, 



£/f (A) ^ $(A) - - , L?(A)-0, 



We recall that fa = Y^keVi r<R(i k) if - - Performing the inverse change of the variables U —> U, 
we get that, for an arbitrary vertex i, 

/•oo 

U?{\) = / e- xt E(J>(X?),X? G L[ iOO |X « = i)dt - *(A) ■ pf (17) 

JO 

where j5f ^ is given by equality (5). For an arbitrary point x G S, 

G L[ iOO |X ™ = x) = E(<j>{X?),T > t\X% = x)l x&LU + 

N ,t 



+ J2 [ E(4>{xt s ),xt s g lf |xj = i)P( T g ds,x? = i), 

i=l 70 



(18) 



where r is the moment when the process X n hits one of the vertices. For x lying on one of 

the rays, the distribution of r coincides with the distribution of the same moment for a Walsh's 

Brownian motion (or a Wiener process on a half-line). On the other hand, it is easy to verify 
p 

that r — >0 for x G So (So = S\S is the joint of all the edges contracting into a point), and this 
convergence is uniform in x G So- Then, for an arbitrary continuous function ip on the limiting 
bundle of rays S, relations (17) and (18) yield 

/•oo 

/ e- xt E(^(X?))\XZ = x)dt^R x i>(0), xGSo, (19) 
J 

e- xt Ety(X?))\XS = x)dt^R x 1>(x), xGS, (19') 



/ 

Jo 
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where R\ip(x) = J °°e~ A * E(ip(X t ))\Xo = x) is the resolvent of the Walsh's Brownian motion X, 
whose asymmetry parameters are given by equality (5). Moreover, the convergence in (19) and 
(19') is uniform on So an d every bounded subset S, respectively. 

We now can present, at last, a martingale characterization of the limiting points of the 
sequence {X n }. Let X be one of such points which is the limit of the subsequence {X Hk }. 
By A, we denote the infinitesimal operator of the above-mentioned Walsh's Brownian motion 
X, whose domain has been described in Proposition 2. For arbitrary A > 0, <fi G Da, to > 
0,ii, . . . ,t m C [0, to],G G Cb(M. m ), by using the Markov property of the processes X Uk and 
relations (19) and (19') for ip = <j) and ip = Acft, we get the equality 

/•oo 

E / e-W-^G(X tl , . . . ,X tm )[\<P(X t ) - A<P{X t )]dt = 

J to 

/•oo 

= lim E e-^-^G^*,...,^)^^)-^^)]^ 



= EG(X tl ,. . .,X tm )[\R^{X t0 ) - R x A<t>(X t0 )} = EG(X tl ,. . . , X tm )cj>{X to ) . 
We now repeat the arguments given in [5], Section 2: relation (20) yields 



(20) 



Xe- xt EG(X tl ,...,X tr: 



to 



cj>{X t )-<t>{X t0 ) 



A<j)(X s )ds 



to 



dt = 0. 



Since a continuous function is uniquely determined by its Laplace transform, we have that 
EG(X tl , • • • , X tm ) \4>{Xt) — 4>{X to ) — f t A<f>(X s ) ds] = 0. Since G and t%, . . .t m are arbitrary, this 
implies that the process 4>{Xt) — f Ac/)(X S ) ds is a martingale. That is, X is a solution of the 
martingale problem presented in Proposition 2. In addition, the distribution of Xq is known 
and coincides with /i. Then, by virtue of Theorem 1, the distribution of X coincides with the 
distribution of a Walsh's Brownian motion with the initial distribution /i. Because we took, as 
X, an arbitrary limiting point of the sequence of distributions {X n }, this means that the whole 
sequence converges to the indicated distribution. Theorem 3 is proved. 



3. Examples 



In this section, we present two examples of the application of Theorem 3. First, consider the 
example, in which the limiting graph 9 consists of one vertex and two half-lines. In this case, 
the limiting process is the well known skew Brownian motion (see j9], j!8j). 

Example 1. Let S consist of two vertices {1,2} and four edges Li,2i ^2,1 , £1,001 and £2,00- 
Let the lengths of the edges L\p and £2,1 be Z™ 2 = n ~ l G and Z 2 i = n~ 1 n, where 8,k > are 
some constants. We suppose that the drift coefficient is zero and the diffusion coefficient is unity 
on all the edges. We write the asymmetry parameters for the vertices {1,2} in the form 

Ploo= l Af, Ph= l ~^f, Ploo= l Af, Ph= 1 ^ (f-*€(-l,l), (21) 

2 = 1,2. Then, by virtue of Theorem 3 (the calculations are simple and omitted), the limiting 
process possesses the asymmetry parameters p± = , where 

1 + q = 9 _ 1 + gi _ l_-q2 (22) 

1 - q k 1 - qi l + ffc' 
Consider two special cases. First, if Z™ 2 = Z 2 i> we have, in essence, the graph with one non- 
oriented edge contracting into a point or, from the other viewpoint, the diffusion process with 
two singular points ("membranes") contracting into a single point. We can rewrite now formula 
(22) as q = th(Ai — A2), Aj = arcthg^, which reproduces the result mentioned in Introduction. 
The difference in signs is caused by the fact that the notation of the asymmetry parameters used 
in [B] differs from (21). 

Another interesting special case arises if qf = 0. In this situation, the prelimiting process is 
a mixture of two Brownian motions on the half-lines (— oo,n _1 0], [— n^ 1 K, +00), for which the 
switch from one motion to another one happens at the moments of hitting of the point n~ 1 8 
(for the first motion) or —n^ 1 K (for the second one). After the switch, a new motion starts from 
the point 0. Such a process is naturally interpreted as the Brownian motion with buffer zones 
on a straight line. 



It is known that the asymmetry of a skew Brown motion at zero can be interpreted in different 
ways: as the result of "a gambling of excursions" with unequal probabilities ([9]), the presence 
of a singular drift coefficient a = qSo ([IE]); or, what is close to the previous, the presence of 
a semipermeable membrane at the point (|12j). Theorem 3 allows us to propose one more, 
apparently completely new interpretation of the skew Brownian motion. Let us consider the 
ordinary Brownian motion (without any singularity!) and introduce the buffer zones in the 
neighborhood of zero with a fixed ratio of the lengths of the zones •& = —. Then, "from the 
macroscopic viewpoint" (i.e., when the size of the zones tends to zero), this motion has the form 
of a skew Brownian motion with the parameter q = = jjr% (see Fig. 1). 



i =- 

1,2 



« M ,« _ K_ 



9-k 

e+K 



Fig. 1. 

Example 2. Let S be a bundle consisting of m half-lines, and let the diffusion process X 
on S have the coefficients of drift and diffusion a r and a r ,r = l,...,m and the asymmetry 
parameters p r ,r = l,...,m. Additionally, we assume that the functions a r ,a r ,r = l,...,m 
satisfy the local Lipschitz condition on R + , so that, from the start from some point on the 
edge till the moment of hitting of the vertex, the process X can be represented as the strong 
solution of an SDE driven by a one-dimensional Wiener process. At the same time, if at least 
three numbers from the collection {p r } are nonzero, then the vertex of the graph S is a triple 
point, and the diffusion X contains the Tsirelson's singularity and, in particular, cannot be 
presented as the strong solution to any system of SDEs driven by a (multidimensional or even 
infinitedimensional) Wiener process. 

We now construct an approximation of this diffusion (in the sense of the convergence in 
distribution) by the sequence of diffusions not containing such a singularity. Consider the graph 
S consisting of m vertices enumerated by the numbers 1, . . . , m and 2m edges. One edge-" ray" 
and one edge-" segment" leave each vertex. The edge-" segment" joins the i-th vertex with the 
i + 1-th one (to shorten the notation, we use the agreement m+ 1 = 1). We put 



i,oo 



a 



a 



a*(0), 



(23) 



Pi,oo — P?,i+l — 2' ^M+l — ~P L -> (24) 

i = 1, . . . , m (we drop the superscript r, because each two vertices are connected by at most one 
edge, and exactly one edge-" ray" leaves each vertex). Let us consider a sequence of diffusion 
processes {^ n } on S with parameters (23), (24) and apply Theorem 3 to it. The limiting graph 
has the single vertex, therefore, we drop the superscript i in the notations given below. We have 



CK*,i+l = 1 and atij = for j ^ i + 1. Therefore, 

fl, j = i + l, 1 

A = 1, Ai j = < , 7Ti = — , 1 = 1... ,171. 

Then, as far as ^2 r p r = 1, the normalizing constant P in formula (5) is P = -jr. Hence, 
the limiting asymmetry parameters are equal to p , ... ,p m . The coefficients a^- , o"™^ , satisfy, 
obviously, the conditions of Theorem 2. Hence, the projections X n of the processes X n on S 
converge in distribution in C(M + , S) to the process X (see Fig. 2). 





Fig. 2. 

By construction, the process X n has no triple points. Let us show that this process can 
be presented as a mixture (in the same sense as that in Example 1) of the collection of strong 
solutions of one-dimensional SDEs. This means that X n can be obtained from a one-dimensional 
Wiener process "in real time". 

Consider m diffusion processes Y™, . . . , Y^, and let the i-th process be defined on a half-line 
L™ = [—^p l ,+oc) and have the coefficients 

ai(x) = a i (0)l [ _i p%0] (x) + a*(z)l+(x), a^x) = ^(0)I h i piQ] (x) + a^x)!^). 

By virtue of the Lipschitz property of the coefficients, each of the processes can be represented 
as the strong solution of a one-dimensional SDE prior to its exit on the boundary of a half-line. 
The process X n can be represented now as follows: at every time moment, the point X n is 
located on one of the half-lines and moves along it by the law defined by the process Y- n 
till the moment of its exit on the boundary of a half-line. Then the point passes to the point 
with the coordinate 0, located on the half-line L™ +1 , and proceed moving along via the law 
of Yj+i, etc. The presented interpretation implies that the filtration generated by the process 
X n can be obtained by a morphism from the filtration generated by a one-dimensional Wiener 
process. That is, the process X n does not contain Tsirelson's singularity. 

Let us summarize: "untwisting" a vertex with multiplicity > 3, i.e., representing a vertex 
with multiplicity > 3 as a result of the "contraction into a point" of a knot, all the vertices of 
which have multiplicity < 2, we have obtained the approximation of a process containing the 
Tsirelson's singularity by processes without such a singularity. It is clear that one can apply the 
same trick to graphs of arbitrary configuration, "untwisting" each vertex with multiplicity > 3. 
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